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ABSTRACT 

We model the optical to X-ray continuum spectrum of the narrow-line Seyfert 1 galaxy 
RE J1034+396. We show that the flat optical spectrum is consistent with emission 
from an irradiated accretion disk. The X-ray emission can be modelled with a disk 
blackbody and a Comptonised component. The temperature at the inner edge of the 
disk T; n = (0.12 ± 0.02) keV. Using this constraint, we show that the flat optical 
spectrum is consistent with emission from the irradiatively heated outer part of the 
accretion disk. We constrain the outer radius of the optically thick disk (i?^5x 10 16 
cm) and the inner radius of the irradiation-dominated region (fl ^ 5 x 10 12 cm). Our 
optical and X-ray spectral fits indicate a mass 0.6 ^ M ^ 3 x 1O 6 M , and do not rule 
out a low (i.e. face-on) inclination angle for the system. 

Key words: accretion, accretion discs — galaxies: active — galaxies: nuclei — galax- 
ies: Seyfert — galaxies: individual: RE J1034+396 — X-rays: galaxies 
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1 INTRODUCTION 

The optical to X-ray continuum of RE J1034+396 is highly 
unusual for a Seyfert 1 galaxy. The optical/UV continua of 
most AGN rise towards the blue with a slope a w 0.4 (where 
q is the spectral index, defined such that F v oc v~ a ), while 
the soft X-ray spectrum falls towards high energies with 
spectral index a w 2 (e.g. Laor et al. 1997). The spectrum 
seems to peak in the unobservable EUV and this continu- 
ous, optical to soft X-ray feature, known as the "big blue 
bump" , is believed to represent the emission from a geomet- 
rically thin, optically thick accretion disk. 

In RE J1034+396, however, the optical/UV continuum 
is flat (q ~ 1), with no sign of the big blue bump down 
to Lya. At ~ 0.1 keV, the soft X-ray spectrum is very 
strong above the extrapolated level of the optical/UV con- 
tinuum, peaking at ~ 0.3 keV then falls steeply towards 
higher energies (Puchnarewicz et al. 1998). Puchnarewicz 
et al. (2001) suggested that the accretion disk component 
in RE J1034+396 is one of the hottest observed in AGN, 
so that its high-energy turnover is measurable in the soft 
X-rays rather than the UV band. They inferred a high ac- 
cretion rate (L ~ 0.3-0.7 Lecm), a small black hole (BH) 
mass M m 1O 6 M0 and an inclination angle i ~ 60-70°. 

The flat optical spectrum was modelled 
by Puchnarewicz et al. (2001) with an underlying power- 
law component extending from the optical to hard X-rays, 
in addition to the optically thick disk spectrum. However, 
they could find no satisfactory physical explanation for this 
feature. Breeveld & Puchnarewicz (1998) detected no sig- 



nificant polarization in the optical, making a synchrotron 
origin unlikely. 

In this Letter, we suggest a different interpretation of 
the continuum spectrum: we show that the flat optical spec- 
trum is consistent with the continuum emission of an irra- 
diated disk, with no need to invoke an additional power-law 
component. We also revise the estimate for the central mass, 
and note that low inclination viewing angles cannot be ruled 
out. 



2 OBSERVATIONS 

Optical spectra of RE J1034+396 were obtained on 1996 
March 24 with the ISIS spectrograph mounted on the 
William Herschel Telescope at La Palma (in service observ- 
ing mode). The spectrum, which covered the range 3800- 
7000 A, is clearly contaminated by light from the host 
galaxy, particularly in the red. In order to remove the ex- 
tended galactic emission, the nuclear spectrum was decon- 
volved from the total spectrum using the technique de- 
scribed in Puchnarewicz et al. (2001). 

In the UV band, RE J1034+396 was observed by HST 
on 1997 January 31 using three gratings (G130L, G190L 
and G270L) covering the range 1100-3300 A. Full details 
of the observations and data reduction were presented in 
Puchnarewicz et al. (1998). 

Beppo-SAX observed the system on 1997 April 18. Use- 
ful data were obtained with the Low Energy Concentrator 
Spectrometer (LECS; 0.1-10 keV) and the Medium Energy 
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Table 1. XSPEC fitting parameters to the Beppo-SAX data 



RE J 1034+393 BeppoSAX (LECS+MECS) 



wabs*(diskbb+pl) 




0.5 1 2 

channel energy (keV) 

REJ1034+393 ASCA (GIS2+3) 




0.5 1 2 5 

channel energy (keV) 

Figure 1. Top panel: disk blackbody plus power law fit to the 
0.15-10 keV Beppo-SAX spectrum. The fit parameters are given 
in Table 1. Bottom panel: the same disk blackbody plus power 
law model (parameters fixed as in Table 1) applied to the 0.4-8 
keV ASCA spectrum. 



Concentrator Spectrometer (MECS; 1.3-10 keV); all three 
MECS units were used. The source was not detected with 
the high-energy instruments. The total observing time was 
21.7 ks for LECS and 43.2 ks for MECS. For more details 
on the observations and on the techniques used for source 
extraction and background subtraction, see Puchnarewicz et 
al. (2001). 



3 X-RAY SPECTRAL ANALYSIS 

In Puchnarewicz et al. (2001), LECS and MECS spectra 
were fitted with a number of models (see their Table 4) in- 
cluding power-law and single-temperature blackbody com- 
ponents. The best fit to the soft spectrum was obtained with 
two blackbody components and a power law, with a column 
density fixed at the Galactic value (uh = 1.5 x 10 20 cm -2 ). 

In this work, we want to test our hypothesis that the 
optical-to- X-ray continuum spectrum is dominated by an 



Parameter 


Value 


wabs*(diskbb - 


h power law): \l = 1-00 (121 dof) 



kh (column density) 


1.5 X 10 20 cm~ 2 


Ti n (T at inner disk edge) 


0.116t° :S£ keV 


.fsfdbb (diskbb normalisation) 


qo'7q + 1368 

J -1068 


r (photon index) 


2 5Q+ - 16 


K p i (power-law normalisation) 


7mt 1 1 1 t x io- 4 



wabs*(diskbb + comptt): xl = 1-05 (119 dof) 



»H (column density) 


1.5 X 10 20 cm" 2 


T; n (T at inner disk edge) 


0.1161°;°°? keV 


K dbb (diskbb normalisation) 




To (input soft photon T) 


o-os±S:Si keV 


T c (T of hot corona) 


14.8 keV 


t (optical depth of hot corona) 


1.307 


N c (normalisation) 


3.54 x 10~ 3 



optically thick accretion disk. This would be analogous to 
Galactic BH candidates in the soft spectral state. Therefore, 
we fitted the Beppo-SAX spectra with a multi-temperature 
disk blackbody plus a power law, using XSPEC version 11.0.1 
(Arnaud 1996). This phenomenological model can be taken 
as an approximation of the spectrum produced by an ac- 
cretion disk plus a hot Comptonising region. The disk is 
responsible for the soft emission; part of the soft photons 
are inverse-Compton scattered to the hard X-ray band as 
they cross the hot corona. We also fitted the data with the 
more physical Comptonisation model of Titarchuk (1994) 
(comptt in XSPEC). 

The best-fit parameters for the two models are listed 
in Tabic 1, and the spectrum fitted with disk blackbody 
plus power law is shown in Figure 1 (top). We fitted the 
LECS and MECS data simultaneously (using 0.13-10 keV 
for LECS, and 0.9-10 keV for MECS), and we applied a con- 
stant factor of 0.9 in the relative normalisation of the two in- 
struments (Fiore et al. 1998). For both the models presented 
here, the best-fit intrinsic column density converged to 0; 
we have therefore fixed the column density to the Galac- 
tic value. For the disk plus power-law model, the power-law 
spectral index a — 1.6 ±0.2 (photon index F = 2.6) is softer 
than for typical AGN (as previously suggested by Pounds, 
Done & Osborne 1995). 

Puchnarewicz et al. (2001) claimed that there is firm ev- 
idence of a hardening in the spectrum at energies £ 7 keV; 
they suggested that a harder power-law component domi- 
nates at higher energies, and that the same power-law com- 
ponent may explain the optical spectrum if it extends to 
those lower energies. A hint of an upturn in the spectrum 
for E £j 7 keV is indeed visible in Figure 1. However, given 
the low signal-to-noise ratio in the higher energy channels, 
and the fact that the upturn is evident only in the LECS 
data (less reliable at those energies), we prefer not to spec- 
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ulate on the physical reasons for the possible existence of a 
hard component. 

We also used the best-fit parameters inferred from the 
Beppo-SAX data (Table 1) to re-fit the ASCA/GIS data 
from 1994 November, previously discussed in Pounds et al. 
(1997) and Puchnarewicz et al. (2001). The result is shown 
in Figure 1 (bottom); we obtain xl = 0.86 (117 d.o.f.) for the 
disk blackbody plus power-law model. Again, although there 
may be marginal evidence for a hardening of the spectrum 
at E 5 keV, the signal-to-noise ratio is too low to allow a 
meaningful interpretation. 

A significant result from both models is the strong con- 
straint on the temperature at the inner radius of the accre- 
tion disk. We obtain kT in = O.U6±^ keV. This is also in 
agreement with earlier ROSAT (Puchnarewicz et al. 1998) 
and ASCA (Pounds et al. 1997; Puchnarewicz et al. 2001) 
observations. We shall use this firm result as we model the 
disk spectrum in the optical band. 



4 CONTINUUM SPECTRUM OF AN 
IRRADIATED DISK 

The effective temperature of an accretion disk photosphere 
at each radius is determined by the thermal energy gener- 
ated by local viscous dissipation near the disk mid-plane 
(viscous heating) , as well as by the energy intercepted from 
the central X-ray source, thermalised and re-emitted (irra- 
diative heating). 

In the absence of irradiation, the effective temperature 
of a geometrically thin, optically thick disk is (Shakura & 
Sunyaev 1973; see also, e.g., Frank, King & Raine 1992) 



TvOR) 



/ 3GMM\ 

\ 8<77r7? 3 J 



1/4 r 



1/2' 



1/4 



(1) 



where R c is a critical radius (< Ri n ) beyond which we as- 
sume no further viscous dissipation in the accretion flow, 
and M is the mass accretion rate at radius R. 

The contribution of the irradiation to the effective tem- 
perature is (e.g., King 1998): 



tL{R) 



r]M c c 2 (l- (3) (H\ n (AYnH 



4naR 2 



/H\ n / ding _ A 
\r) \d\nR )' 



(2) 



where r\ is the efficiency of conversion of accreted mass into 
energy, M c is the central mass accretion rate (M c < M(R) if 
some mass is lost in a disk wind) , /3 is the albedo of the disk 
surface, H is the disk scale-height at radius R. The index n 
(1 < n < 2) depends on the geometry of accretion (Vrtilek 
et al. 1990; van Paradijs 1996). 

The last two factors on the right hand side of equation 
(2) depend on the disk thickness and the amount of disk 
flaring (a more concave disk intercepts a larger fraction of 
radiation). They vary only slowly with R; following King 
(1998), we approximate them with a constant factor 7 = 
(H/R) n (dlnH/d\nR - 1), with 0.01 £ 7 £ 0.05 for a thin 
disk. Hence, T irr ~ R~ 1/2 while T v ~ i?~ 3/4 ; this implies 
that irradiative heating dominates at large radii, and viscous 
heating at small radii. 

Henceforth, we shall also assume that M c — M(R), and 
re-write the mass accretion rate in terms of the emitted lu- 
minosity: M(R) = (L/L Edd ) (L Edd /ric 2 ) = l Edd (L E dd/V 2 ), 
with L Ed d = 1.3 x 10 38 (M/Mq) erg s -1 . Finally, given the 



Table 2. Best-fit parameters for the irradiated disk model [sec 
Eqs. (3), (7), (8)] 



Parameter 



Value 



1.5 

7.5 x 10 5 K 
43.5 

0.56 x 10 12 cm 
5.0 x 10 16 cm 
0.235 



characteristic dimensions of RE J1034+396, we shall define 
n 2 = i?/(10 12 cm) and m 6 = M/(1O 6 M ). 

From equations (1) and (2), the total effective temper- 
ature is given by: 



T cS (R) = [T*(R)+tUR)] ar 
with the parameters a and I given by: 



-3/4 



1/4 



a ~ 6.7 x 10" 



0.2 



1/4 



,1/4 1/2 

f Edd m 6 



K 



, (3) 



(4) 



Hence, we have 



T in = T cB (R in ) ~ a ( 



Rin 



10 12 cm 



-3/4 



-3/4 



(6) 



The specific flux at frequency v detected by an observer 
at distance d, with the line of sight forming an angle i to 
the normal to the disk plane, is (Mitsuda et al. 1984): 



F v 



47rfe(cos i)v 
c 2 d 2 (l + «) 2 



RdR 



exp (hv/kT eS ) 



(7) 



where the redshift 1 + z = 1.042 for RE J1034+396. This 
corresponds to a proper distance d = 3.8 x 10 26 /i -1 cm. Here 
we have defined, as usual, h = Ho/(WO km s~ Mpc~ ). 

Equation (7) is sometimes modified with the introduc- 
tion of a spectral hardening factor / (Ebisuzaki, Hanawa 
& Sugimoto 1984; Shimura & Takahara 1995), to take into 
account Compton scattering in the disk. We can write the 
"diluted disk blackbody" specific flux as: 



F* h = r 4 F v {fT eS ) 



(8) 



with 1 < / < 1.7 (Shimura & Takahara 1995) 

Our goal is to check if a disk blackbody or a diluted 
disk blackbody spectrum can explain the observed opti- 
cal/UV/soft X-ray continuum. We have compared the ob- 
served fluxes vF v with the values predicted by equations (7) 
and (8), for various values of the parameters 
/ and i. Using the results of our X-ray spectral fitting, we 
constrain Ti n to vary in the range 0.10 < kT ln < 0.15 keV. 
We shall discuss our main results in the next section. 
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5 DISCUSSION OF OUR RESULTS 

We show in Figure 2 that the continuum spectrum at ener- 
gies ^ 0.3 keV can be well fitted by a diluted disk blackbody 
spectrum (thick solid line) with / = 1.5, and best-fit param- 
eters listed in Table 2. For E £ 0.5 keV, the thick solid line 
plotted in Figure 2 represents the best fit to the Beppo-SAX 
X-ray data, discussed in Section 3. 

The model parameters in Table 2 imply kTi n = 0.10 
keV. For h = 0.5, it is i = 20°, while for h = 0.75, i = 60°. 
Therefore, unlike what was claimed in Puchnarewicz et al. 
(2001), we cannot rule out that the system is seen face-on 
(i 45°). This would explain why the molecular torus does 
not block out the X-ray emission from the central source. It 
has been suggested (Puchnarewicz et al. 1992; Boiler et al. 
1997) that narrow-line Seyfert 1 galaxies are generally seen 
face-on. 

If we assume a simple disk blackbody spectrum instead 
(/ = 1), it is clear from equations (7) and (8) that the same 
fit can be obtained with a = 1.1 X 10 6 K, h 2 cosi = 0.04 
(i — 80° for h — 0.5), with the other parameters fixed as 
above. This corresponds to kT[ n = 0.15 keV. In fact, we 
can obtain good fits for 0.10 fcTin 0.15 keV by varying / 
between 1.5 and 1. Given the uncertainty in the distance and 
in the hardening factor, we cannot constrain the inclination 
angle more strongly. 

The parameter I provides the size of the disk region 
dominated by irradiative heating. In our best fit, the disk 
is irradiatively dominated for R ^ 80Ri n . Small changes in 
the parameter I have large effects on the predicted optical 
spectrum. In Figure 2, we also plot the spectrum obtained 
with I = 60 (marked with "b"; irradiation dominating for 
R ^ 110 -Rin) and I = 30 (marked with "c"; irradiation dom- 
inating for R ^ 50i?in). Using equation (5), we see that the 
best-fit value of I is consistent with M ~ 1O 6 M0, as ex- 
pected. 

We can obtain a stronger constraint on the central mass 
by inserting our best-fit value of a into equation (4): 

M ~ 1.25 xl0 6 ; Ec ^ 2 (^) M . (9) 

For the high accretion rates observed in this system 
(0.3 ^ tedd ~ 1, see e.g. Puchnarewicz et al. 2001), and for 
0.05 ^ r\ £ 0.4, we can constrain the mass of the central ob- 
ject to be 0.6 x 10 6 < M < 3 x 10 6 M Q . 

A simple physical interpretation of Ri n is that it corre- 
sponds to the radius of the innermost stable circular orbit 
around the central BH. However, the standard disk model 
is probably not applicable near the inner disk radius, and in 
any case it has to be modified with relativistic terms (e.g., 
Novikov & Thorne 1973; Riffert & Herold 1995). Therefore, 
our spectral fit near the peak of the blackbody component 
(at ~ 0.1-0.5 keV) can only be a general approximation. 
A more exact modelling is left to further work, and will 
require more precise observational data in the soft X-ray 
band. Nevertheless, we can check that our best-fit value of 
J?in = 0.56 x 10 12 cm is consistent with the range of values 
for M derived before. If M < 0.65 x 10 6 M Q , R in is consistent 
with the innermost stable circular orbit of a Schwarzschild 
BH, while if M > 1.5 x 1O 6 M we require a Kerr BH with 
angular parameter a > 0.85. 

The outer disk radius is also constrained by our fit: we 




14 15 16 17 18 

og v (Hz) 



Figure 2. Broadband spectral fit to the WHT, HST, and Beppo- 
SAX observed fluxes (E < 7 keV). (The line-of-sight Galactic ab- 
sorption has been removed.) The thick solid line for log v < 16.7 
represents the low-energy part of an irradiated disk blackbody 
spectrum, continued as a dash-dotted line for logy > 16.7. The 
thick solid line for logy > 16.9 represents the X-ray spectral fit, 
physically interpreted as a disk blackbody modified by Compton- 
isation. The temperature at the inner edge of the disk fcT; n = 0.1 
keV. The meaning of the curves labelled "a" , "b" and "c" is dis- 
cussed in the text. 

require Rout ^5x 10 16 cm « 3 x 10 5 GM/c 2 . If we keep 
all other parameters in Table 2 fixed, but assume a disk 
truncated at R ou t = 10 16 cm, we cannot reproduce the flat 
optical spectrum (curve marked "a" in Figure 2). 



6 FURTHER WORK: DISK WINDS AND LINE 
EMISSION 

A large accretion disk heated by soft X-rays intercepted from 
the central source is likely to have a temperature-inversion 
layer at its surface (Wu et al. 2001). Therefore, we expect 
strong emission-line formation at large radii. In particular, 
in the case of RE J1034+396, the largest contribution to low- 
ionisation lines such as Ha and H/3 would come from radii 
^ 10 16 cm, where the Keplerian rotational velocities are 
(GM/R) 1/2 < 1000 km s -1 . If the Balmer lines are indeed 
produced near the disk surface, we expect full widths at half 
maximum ~ 2(sini)(GM/R) 1/2 < 2000 km s" 1 . This is in 
agreement with the observed values of 1500 and 1800 km 
s _1 for H/3 and Ha respectively (Puchnarewicz et al. 1995). 

Moreover, if the X-ray luminosity from the central 
source is 0.1 L&dd, we may expect the formation of a 
radiatively-driven accretion disk wind. Accurate modelling 
of the emission-line profiles may help us ascertain if higher- 
ionisation lines are formed in the wind, and determine the 
ionisation parameter of the emission regions. The narrow, 
single-peaked UV lines seen by Puchnarewicz et al. (1998) 
appear qualitatively consistent with the disk-wind model of 
Murray & Chiang (1997). The presence of broad and nar- 
row components in the UV lines may be explained with a 
broad component emitted near the irradiated disk surface at 
smaller radii, and a narrow component emitted in the pho- 
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toionised wind. A full investigation of this issue is beyond Wu, K., Soria, R., Hunstead, R. W., Johnston, H. M., 2001, MN- 

the scope of this Letter and is left to further work. RAS, 320, 177 



7 CONCLUSIONS 

We have fitted the optical to X-ray continuum spectrum 
of the narrow-line Seyfert 1 RE J1034+396 with an irra- 
diated disk blackbody component (temperature kT[ n m 0.1 
keV at the inner disk radius), and a power- law- like Comp- 
tonised component at higher energies. The disk component 
is much stronger and found at higher energies than in the 
average AGN; the spectral index of the power-law compo- 
nent is a ~ 1.6 (steeper than the average value for an AGN); 
the observed luminosity is close to the Eddington limit. We 
infer a central mass of 0.6 x 10 6 < M < 3 x 10 6 M Q . 

The observed flat optical/UV spectrum is explained by 
our model as emission from the radiatively-heated outer part 
of the disk, at R > 80 R in , where T ~ R~ 1/2 . Strong emis- 
sion lines are likely to be formed near the irradiated sur- 
face of the outer disk, and in a radiatively-driven wind. Our 
model requires the disk to extend to R out w 5 x 10 16 cm. 
The observed flux and inner-disk temperature do not rule 
out the hypothesis of a face-on system, which would be in 
agreement with the low line-of-sight absorption seen in the 
X-ray spectra. 
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